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Abstract 


Models  based  on  empirical  fluid  flow  and  heat  transfer  data  combined  with  basic  principle 
analysis  provide  data  for  estimating  the  cost  of  implementing  a  closed-loop  radiant  cooling 
system  for  ships  in  dry  docks.  Depending  on  the  material  used,  pipe  lengths  from  300  m  to 
3000  m  with  diameters  ranging  from  10  in  to  24  in  all  cooled  a  design  load  of  8  F°  = 

4. 44  K  above  ambient  temperature  when  that  temperature  was  assumed  to  be  300  K.  From 
the  cost  analysis  performed,  $16  million  would  be  the  maximum  initial  investment  for  the 
system  that  moves  26  million  gallons  per  day  (~  18000  gpm).  Although  this  analysis  has 
not  accounted  for  redundancy,  pumps,  and  other  one-time  installation  costs,  other  solutions 
for  this  problem  have  proposed  budgets  of  $40  million,  making  this  solution  cost- 
competitive. 
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Foreword  [1] 

Most  Navy  vessels  in  dry  dock  only  partially  deactivate.  Many  of  their  components 
continue  to  operate  normally  and  require  cooling  water.  Currently,  this  water  is  simply 
pumped  from  the  harbor,  through  the  ship's  systems,  and  discharged.  The  effluent  cooling 
water  contains  a  higher  concentration  of  copper  and  other  metals  and  is  at  a  higher 
temperature  compared  to  ambient  conditions;  this  is  not  a  concern  at  sea,  but  at  dock 
additional  environmental  regulations  apply.  These  regulations  will  become  stricter  in  the 
short  to  medium  term,  and  this  will  require  a  change  in  current  practice.  A  closed-loop 
water  cooling  system  would  eliminate  the  discharge  of  any  cooling  water  back  into  the 
harbor  and  thereby  satisfy  regulatory  requirements. 

A  closed-loop  cooling  system  solves  this  problem  by  continuously  recirculating  the 
same  cooling  water  through  a  dry-docked  ship.  Since  no  water  is  discharged  into  the  harbor, 
environmental  regulations  are  satisfied.  The  challenges  with  closed-loop  cooling  systems 
based  on  traditional  engineering  solutions  like  industrial -scale  refrigerator  or  chilling  are 
scaling  the  system  to  accommodate  the  volume  of  cooling  water  required  (projected  to  be 
up  to  18  million  gallons  per  day  for  Ford-class  carriers)  without  making  the  physical 
footprint  of  the  system  prohibitively  large,  and  the  ongoing  operating  costs,  in  both 
electricity  and  maintenance,  required  for  such  a  system. 

We  propose  to  address  both  of  these  challenges  through  a  novel  reversal  of 
commercial  technology  widely  used  in  home  and  office  buildings:  radiant  floor  heating. 


Summary 

The  objective  of  this  project  is  to  determine  if  radiant  cooling  technology  is  cost- 
competitive  and  feasible  enough  for  deployment  at  a  Navy  shipyard.  Selection  of  various 
diameters  and  flow  rates  determined  the  required  lengths  of  pipe  to  cool  a  8  F°  load  and  the 
respective  pressure  losses.  These  lengths  and  pressure  losses  specified  the  system  costs  at 
various  flow  rates  for  comparison.  When  plotted  for  various  diameters  and  numbers  of 
modules,  the  costs  associated  with  making  the  system  to  satisfy  a  26  million  gallon  per  day 
flow  depict  a  tradeoff  between  initial  and  recurrent  costs.  Twelve  possible  installation  sites 
at  Puget  Sound  Naval  Shipyard  can  house  the  proposed  system.  Lastly,  the  results  show  that 
the  system  will  cost  $16  million  in  material  costs.  Although  this  analysis  has  not  accounted 
for  redundancy,  pumps,  and  other  one-time  installation  costs,  other  solutions  for  this 
problem  have  proposed  budgets  of  $40  million.  Basic  cost  estimates  for  the  unaccounted- 
for  work  suggest  that  the  radiant  cooling  solution  is  cost-competitive. 
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Introduction 

The  subjects  addressed  here  include  cooling,  fluid  dynamics,  and  Pareto  analysis. 

Cooling  is  a  subset  of  the  discipline  of  heat  transfer  in  which  energy  is  extracted 
from  a  system  as  heat.  In  this  subset,  two  modes  of  cooling  will  be  addressed:  conduction 
and  convection.  Conduction  results  from  physical  contact  of  media.  Convection  arises  from 
the  bulk  movement  of  a  medium.  A  heat  transfer  coefficient,  defined  as  the  amount  of  heat 
flux  experienced  per  a  known  temperature  difference,  captures  the  effectiveness  of  a 
convective  heat  transfer  process.  For  the  case  of  this  problem,  convection  occurs  from  a 
fluid  flow  in  a  pipe  and  conduction  occurs  in  the  soil  surrounding  the  pipe. 

Despite  having  many  analytic  models,  convection  is  mainly  an  empirical  or  numeric 
science.  This  is  especially  true  for  turbulent  flow  where  the  bulk  motion  of  fluid  particles  is 
chaotic  and  involves  a  lot  movement  besides  the  bulk  motion  of  the  fluid.  Since  this 
application  will  consistently  have  convection  from  a  turbulent  fluid  flow,  empirical 
equations  from  past  literature  [6]  will  model  the  convective  heat  transfer.  These 
relationships  give  a  relationship  between  the  heat  transfer  coefficient  and  known  properties 
of  the  flow  being  analyzed.  After  the  heat  transfer  coefficient  is  found,  the  heat  transfer  rate 
of  the  flow  can  be  predicted  within  a  certain  error. 

Fluid  dynamics  becomes  important  any  time  a  flowing  fluid  is  examined.  As  with 
convection,  the  dynamics  and  behavior  of  turbulent  flow  is  mainly  an  empirical  or  numeric 
science.  The  information  of  interest  here  is  head  loss  of  the  flow.  Head  loss  relates  to  the 
pressure  lost  from  a  flowing  flow.  It  can  be  represented  by  a  pressure  or  the  height  of  a 
column  of  the  fluid  since  a  column  of  fluid  exerts  a  pressure  proportional  to  its  density  and 
the  gravitational  acceleration  (i.e.  p  =  pgh).  With  a  known  head  loss,  a  pump  can  be  sized, 
or  chosen,  so  that  it  can  maintain  the  desired  flow. 

A  Pareto  chart  shows  options  that  cannot  have  improvement  in  one  area  without 
causing  a  decline  in  another.  Such  a  chart  reveals  the  existence  of  a  tradeoff  between 
components  of  the  options  being  compared  without  a  way  of  knowing  the  value  of  each 
area  and  how  the  values  relate.  All  shown  points  on  a  Pareto  chart  are  the  most  relatively 
efficient  options  in  the  space  depicted  by  the  axes  used  to  plot  the  chart. 

Notkin  Engineering  previously  did  an  extensive  engineering  report  on  solutions  to 
the  dry  dock  problem.  The  report  addressed  many  solutions  that  did  not  include  this  radiant 
cooling  option.  In  the  conclusion,  adjusting  current  infrastructure  to  meet  stricter 
regulations  and  still  cool  the  ships  solved  the  problem  most  cost-effectively.  Improvements 
cost  an  estimated  $42  million  at  the  least  and  $49  million  at  the  most.  [3] 

The  purpose  of  this  investigation  is  to  determine  the  cost-competitiveness  and 
feasibility  of  implementing  a  system  that  cools  a  fluid  passively  in  the  ground  for  use  in  dry 
docks.  First,  the  resulting  equations  from  analysis  are  presented  with  the  assumptions  used. 
Second,  the  process  for  determining  the  cost  of  the  system  is  discussed.  Third,  the  results  of 
the  equations  are  displayed  graphically  and  cost  comparisons  are  made.  Fourth,  areas  in 
Puget  Sound  Naval  Base  that  provide  enough  space  for  construction  are  identified.  Finally, 
a  recommendation  with  extra  thoughts  will  be  provided  on  what  option  to  use.  The  scope  of 
this  exploration  encompasses  the  cost  associated  with  pipe  and  estimated  power  costs  for 
the  system,  not  the  cost  of  the  emergency  redundancy,  pumps,  excavation,  and  planting  of 
the  system  in  the  ground.  From  this,  managers  and  engineers  working  in  dry  docks  can 
determine  whether  this  system  would  be  a  wise  investment  and  what  pumps  to  use  later. 
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Methodology 

Two  analytic  methods  determined  the  length  required  for  the  necessary  heat  transfer 
of  the  problem  and  the  resulting  pressure  loss  from  this  length  when  bent  into  a  radiator  like 
form.  All  calculations  used  SI  units.  When  English  units  were  the  desired  inputs  or  outputs, 
their  metric  equivalents  were  actually  calculated  and  then  converted.  Appendix  A  holds  the 
detailed  derivations  of  the  equations. 

For  the  heat  transfer  analysis,  assumptions  made  computation  tractable.  First 
assumption:  the  fluid  in  the  pipe  flowed  in  a  steady,  fully-developed,  incompressible 
stream.  As  the  interest  is  in  the  performance  of  a  turbulent  flow  in  steady  state,  steady  flow 
is  a  valid  assumption,  and  the  flow  would  reach  full  development  quickly.  The  fluid  will 
consist  of  a  liquid,  and  most  liquids  compress  negligibly.  Second  assumption:  the  flow 
maintains  a  constant  average  velocity  since  the  pipe  diameter  is  to  remain  constant.  Third 
assumption:  the  system  losses  energy  only  by  heat  transfer.  This  is  sensible  because  the 
system  will  have  no  work  extracting  components.  Fourth  assumption:  all  friction  factors  are 
Darcy  friction  factors  as  it  is  the  most  commonly  used  parameter  in  the  derivations  of  the 
pipe  flow  equations.  Fifth  assumption:  the  characteristic  length  of  the  system  is  the  pipe 
diameter  because  the  equations  used  follow  this.  Lastly,  the  pipe  will  initially  be  highly 
conductive  so  that  the  conduction  resistance  across  the  pipe  could  be  neglected. 

Fig.  1  shows  the  simplified  diagram  for  the  heat  transfer  analysis. 
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Fig.  1  -  Simple  Heat  Transfer  Model  Sketch 

Using  this  figure,  a  simple  energy  balance  found  the  heat  lost  as  a  product  of  the  Nusselt 
number,  surface  temperature  gradient,  and  surface  area.  To  evaluate  the  surface  temperature 
gradient,  the  average  bulk  temperature  of  the  fluid  was  used  and  the  surface  was  held  at 
constant  temperature.  After  some  algebra,  the  length  of  the  pipe  follows  equation  (1): 

^  _  pD^VCpO\-T2) 

_  2(,Nuk)(T1+T2-2Ts)  ^  ' 

where  the  Nusselt  number  is  evaluated  using  an  empirical  relationship.  If  the  pipe  was  not 
conductive,  the  heat  lost  would  be  represented  as  a  function  of  the  same  product  as  before 
but  with  the  Nusselt  number  term  replaced  by  the  reciprocal  of  the  sum  of  the  thermal 
resistances  of  the  convective  transfer  and  conductive  transfer.  This  results  in  equation  (2): 

J  _  F/4+V(F^)2+4FBCDE7T 
_  2CDEn  ' 

A  =  2nDkp;  B  =  Nu  /cln((D  +  2t)/D);  C  =  2nkpNu  k; 

E  =  (T±  +  T2-  2Ts);  F  =  P-^cp(T1  -  T2) 

For  the  pressure  drop  analysis,  two  additional  assumptions  are  made.  First  new 
assumption:  the  system  remains  at  a  constant  elevation  in  the  ground.  This  is  reasonable 
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since  the  system  is  expected  to  be  installed  horizontally.  Second  new  assumption:  the 
system  will  consist  of  a  known  number  of  180°  bends  and  straight  sections  of  pipe  denoted 
by  /?  and  y  respectively. 

The  pressure  drop  analysis  started  with  a  mechanical  energy  balance  and  simplified 
to  the  following  relationship: 

Ap=p(yK  +  (]f±)  |  (3) 

With  both  of  these  relationships,  various  design  values  gave  length  and  pressure  loss 
values  for  various  diameters  (10  in  —  24  in),  pipe  materials  (conductive  and  insulative), 
and  flow  rates.  Conductive  materials  had  thermal  conductivities  above  15  W /mK.  A  flow 
rate  of  interest  was  26  million  gallons  per  day.  The  knowns  of  the  relationships  came  from 
a  working  fluid  of  25%  sodium  chloride  (NaCl)  aqueous  solution,  a  bend  of  R  =  4.5 D  (D  is 
the  diameter  of  the  pipe),  a  smooth  schedule  40  pipe  (< e  =  0),  a  cooling  load  of  8  F°  above 
soil  temperature,  a  soil  temperature  of  300X,  and  the  system  having  5  bends  and  6  straight 
sections.  These  knowns  can  be  updated  in  the  Python  notebook  in  Appendix  B  to  regenerate 
the  analysis  for  the  new  information. 

Costs  were  determined  from  the  results  of  the  two  analyses.  The  total  initial  cost  of 
the  system  was  calculated  as  the  cost  of  the  required  length  of  pipe.  The  ongoing  cost  was 
kept  in  terms  of  the  power  usage  of  the  system  for  easy  conversion  later.  Using  the  lengths 
calculated,  the  cost  of  insulative  piping  came  from  the  typical  price  values  of  PVC  pipe  at 
various  diameters  advertised  in  2015  [4],  Conductive  piping  costs  were  found  for  10  in 
and  12  in  diameters  [7].  The  rest  were  estimated  by  taking  the  insulative  costs  and 
multiplying  by  the  ratio  of  metal  pipe  cost  and  plastic  pipe  at  12  in  diameter.  Appendix  D 
summarizes  the  prices  used.  The  power  required  for  the  flow  was  calculated  using  the  pump 
power  equation  assuming  a  pump  efficiency  ofr]  =  82%  with  the  equation  P  =  VAPr], 

The  above  costs  cover  implementation  of  one  big  system  to  sate  the  needed  flow 
rates.  In  practice,  many  smaller  yet  similar  systems,  called  modules,  work  together  to  cover 
the  desired  flow  rates.  This  increases  redundancy  and  improves  the  pressure  head.  With  the 
division  into  modules,  each  module  will  experience  less  total  flow  rate  affecting  the 
required  lengths  and  power  costs.  Total  costs  for  a  modularized  system  are  the  cost  of  each 
module  times  the  number  of  modules. 

There  already  existed  a  tradeoff  between  power  and  material  costs  but  now  there  is 
also  a  dependence  on  how  many  modules  will  be  used.  A  Pareto  Chart  will  depict  the 
tradeoff  present  when  evaluating  costs  at  the  flow  of  interest  (26  gallons  per  day),  various 
diameters,  and  various  numbers  of  modules.  Once  this  chart  is  created,  various  cost 
estimates  will  help  in  comparing  costs. 

To  determine  where  the  system  could  be  implemented,  the  total  area  was  estimated. 
With  the  knowns  above,  a  rectangle  could  enclose  the  design  according  to  equation  (4): 

A  =  (a~^7rfl)  +  2 R)  (y2R)  (4) 

A  large  length  of  about  1650  m  was  used  to  get  a  rectangle  that  was  big  enough  for  most 
cases.  The  amount  of  area  was  located  throughout  the  Puget  Sound  Naval  Base  in  space 
where  there  was  only  continuous  parking  lot,  grass,  or  gravel  on  the  main  land.  The  area 
was  assumed  to  be  malleable,  so  the  original  rectangle  was  not  always  maintained  but  the 
area  remained  essentially  the  same.  A  different  color  designates  each  separate  viable 
location. 
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Results  and  Discussion 

The  heat  analysis  equations  gave  graphical  outputs  of  Fig.  2  and  3. 

Length  vs.  Volume  Flow  Rate  in  Highly  Conductive  Pipe 
Volume  Flow  Rate  (gpm) 

5000  10000  15000 


Fig.  2  -  Required  Length  vs.  Flow  Rate  -  Conductive  Material 
Length  vs.  Volume  Flow  Rate  in  Insulating  Pipe 
Volume  Flow  Rate  (gpm) 

5000  10000  15000 
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A  length  range  from  300  m  to  2200  m  contains  all  the  observed  possibilities.  All  curves 
share  the  same  trends  that  length  increases  with  flow  rate  and  bigger  pipe  diameters  result 
in  longer  pipes  at  lower  flow  rates. 

The  equations  used  in  developing  the  heat  transfer  analysis  will  carry  on  the  errors 
associated  with  the  values  they  use.  If  the  fluid  property  and  geometry  values  can  be 
considered  as  absolute  constants  (i.e.  no  error  associated  with  them),  the  error  in  the  length 
will  equal  the  error,  in  percent,  of  the  calculated  Nusselt  number.  This  would  give  an  error 
of  less  than  one  percent  for  the  flows  observed  [2],  The  crossing  of  curves  for  the  insulative 
pipes  occurs  due  to  the  subtle  effect  of  the  thickness  variations  for  different  diameters  of 
schedule  40  pipe. 

The  pressure  analysis  equations  gave  graphical  outputs  of  Fig.  Cl  and  C2  located  in 
Appendix  C.  As  a  general  trend,  the  pressure  grows  rapidly  with  increasing  flow  rates, 
following  a  square  function  (i.e.  x2)  to  be  exact.  Also  pressure  increases  with  decreasing 
diameters.  With  the  same  concessions  as  with  the  heat  transfer,  the  pressure  loss  calculation 
gets  its  error  from  the  friction  factor  calculation  [6]  and  the  length.  This  results  in  around  a 
2.25%  error. 

The  cost  analysis  gave  Fig.  C3,  C4,  C5,  and  C6,  located  in  Appendix  C.  The  costs 
for  the  pipe  inherit  the  error  of  1%  from  the  length  calculation  while  the  power  cost 
receives  the  2.25%  error  from  the  pressure  values.  Fig.  C7,  located  in  Appendix  C  was 
created  from  combining  the  values  found  above  evaluated  at  a  flow  rate  of  26  million 
gallons  per  day,  various  diameters,  and  various  module  numbers.  When  Fig.  C7,  had 
decidedly  worse  options  removed,  Fig.  4,  a  Pareto  Chart,  results: 


After  noting  that  this  graph  has  the  errors  from  the  different  costs,  the  most  expensive 
option  initially  is  less  than  $16  million. 
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Twelve  million  dollars  does  not  account  for  emergency  redundancy,  pump  costs, 
and  excavation  costs.  However,  as  a  naive  estimates,  pump  costs  could  total  to  around  $1.6 
million  for  16  pumps  at  $100,000  a  pump.  The  cost  of  digging  is  estimated  to  around  $14 
million  [5],  To  estimate  the  cost  of  redundancy,  10%  will  be  added  to  the  total  so  far  (i.e. 
$3.16  million).  This  projects  the  total  of  the  system  to  cost  around  $35  million  compared 
to  the  $42  million  proposed  by  Notkin  Engineering.  The  cost  of  a  specific  system  is  broken 
down  in  Appendix  E  resulting  in  about  $27  million. 

A  total  area  of  7950  m2  would  contain  the  system.  The  12  spots  found  in  Puget 
Sound  Naval  Base  are  shown  on  a  Google  Earth  image  in  Fig.  5. 


Fig.  5  -  Viable  Area  Locations 


As  Fig.  5  shows,  there  are  many  valid  locations.  However,  all  of  these  are  not  in  the  direct 
vicinity  of  the  dry  docks  to  be  serviced  and  would  require  piping  to  connect  the  areas  with 
the  dry  dock. 

Conclusions 

The  cost  of  the  system  to  cool  an  8  F°  load  in  26  gallons  per  day  with  a  25%  NaCl 
aqueous  solution  in  a  pipe  system  would  cost  around  $16  million.  Although  this  analysis 
has  not  accounted  for  redundancy,  pumps,  and  other  one-time  installation  costs,  other 
solutions  for  this  problem  have  proposed  budgets  of  $40  million.  Basic  cost  estimates  for 
the  unaccounted-for  work  suggest  that  the  radiant  cooling  solution  totals  around  $27  million 
for  a  specific  system.  There  exist  12  areas  where  this  system  can  be  implemented.  Overall 
this  system  is  both  cost-effective  and  feasible. 

Recommendations 

This  analysis  has  shown  that  a  radiant  cooling  system  is  a  cost -competitive  option  to 
address  stricter  environmental  restrictions  and  larger  cooling  demands  at  dry  docks.  This  is 
an  expensive  option  but  possibly  cheaper  than  other  proposed  ones.  Overall,  this  report  has 
shown  that  this  technology  should  at  least  be  considered  for  implementation  at  dry  docks. 
The  cost  estimates  should  be  considered  as  naive  estimates.  Further  research  or  contracting 
can  look  into  true  costs  associated  with  implementing  this  system  in  Puget  Sound  Naval 
Base. 
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Appendix  A  -  In-Depth  Analyses 

Pipe  Length  Analysis  for  Pipe  System 

Description:  A  level  pipe  system  (ps)  is  proposed  as  to  carry  water  in  a  radiant  cooling 
project  design.  As  a  modeling  task,  the  length  required  for  reaching  a  specified  temperature 
in  a  simplified  setup,  i.e.  along  a  straight  pipe  in  the  ground,  is  to  be  analyzed.  The  goal  is 
to  find  a  relationship  between  l,  D,  and  heat  lost. 

Given: 

D  -  diameter  of  the  pipe  system 
v  -  average  speed  of  fluid  in  the  ps 
p  -  the  density  of  working  fluid 
p  -  dynamic  viscosity  of  working  fluid 
Ti  -  inlet  design  temperature 

T2  -  outlet  design  temperature 

k  -  thermal  conductivity  of  working  fluid 
a  -  thermal  diffusivity  of  working  fluid 
Ts  -  temperature  of  the  soil  surrounding  the  pipe 
Determine:  The  length  of  pipe  (0  required  for  given  design  temperatures 
Model:  1)  Steady,  fully  developed,  incompressible  flow 

2)  Constant  average  speed 

3)  Turbulent  flow 

4)  Heat  transfer  is  sole  mode  of  energy  loss 

5)  Darcy  friction  factors 

6)  Highly  conductive  pipe 

Basic  Equations:  Reynolds  Number  Equation 

Re  =  (1) 

Prandtl  Number 

Pr  =  -  (2) 

a 

Nusselt  Number 


nr  hD 

Nu  =  — 

k 


Nusselt  Number  Empirical  Equation  by  Gnielinsk 


Nu  = 


(//8)(fie-103)Pr 


(3) 

(3000  <  Re  <  106 
0.5  <  Pr 


l 


Nu  =  0.023  Re0  8Prr 


l  +  12.7(//8)1/3(Pr2/3-l) 

Nusselt  Number  Empirical  Equation  2 

104  <  Re 
0.5  <  Pr 

n  =  0.33  for  cooling 1 
n  =  0.4  for  heating 
General  Heat  Loss  Equation  (similar  to  Newton’s  Law  of  Cooling) 

Q  =  hAAT  (6) 

Steady  State  Conservation  of  Energy  Equation 


(4) 

(5) 


1  https://en.wikipedia.org/wiki/Heat  transfer  coefficient  referencing  F.  Kreith,  ed.  (2000).  The  CRC 

Handbook  of  Thermal  Engineering.  CRC  Press. 
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Q  +  W  =  mcpAT  (7) 

Analysis: 

First  use  (1)  and  (2)  with  the  given  information  to  get  Re  and  Pr.  From  this  determine 
whether  to  use  (4)  or  (5)  to  calculate  the  Nusselt  Number.  Using  (4)  requires  knowledge  of 
the  Friction  Factor  of  the  flow.  The  Friction  Factor  coefficient  is  a  nondimensional  term 
dependent  on  the  Reynold’s  number  and  Surface  Roughness  coefficient  (-).  With 

assumption  1,  it  can  be  found  using  a  Moody  Plot  or  an  empirical  fitting  equation  like  the 
Colebrook  equation: 


Now  solve  for  the  Nusselt  Number  for  use  in  (3)  to  find  the  heat  transfer  coefficient  of  the 
process: 


Next  (7)  is  simplified  with  assumption  4  and  (6)  and  (8).  Equation  (6)  uses  the 
temperature  of  the  soil  and  the  average  temperature  for  the  temperature  gradient  causing  the 
heat  flow  out  of  the  pipe.  The  use  of  soil  temperature  is  justified  by  assumption  6: 


(77)  A(Tave  -  Ts)  =  mcp(T1  -  T2) 
=>  (^)  nDl(Tave  -  Ts)  =  P-^vcv{Jx  -  T2] 
=>  ( Nu  k)l  (7  +  7  ~TS)  =  P~^vcv(Jx  -  T2) 


=>  ( Nu  k)l(Ti  +  T2-  2 Ts)  =  P-^vcv{Jx  -  T2) 

=>  l  =  pdZ'°pK-V  (io) 

2(vufc)(r1+r2-27’s)  v  ’ 

Equation  (10)  gives  the  length  required  for  the  desired  temperature  difference  for  a 
characterized  flow  in  a  straight  pipe. 

Discussion: 


One  of  the  major  factors  that  made  this  solvable  was  choosing  the  characteristic  length  of 
the  flow  to  be  the  diameter  of  the  pipe  not  the  to-be-solved  for  length  of  the  pipe.  This 
analysis  provides  a  relationship  between  the  length,  diameter,  and  a  desired  heat  loss 
through  design  temperatures.  The  cost  of  computing  this  simple  result  will  be  in  finding  / 
and  calculating  Nu. 

The  assumptions  used  in  this  analysis  are  justified.  The  steady,  fully-developed  flow 
aligns  with  the  desire  for  a  steady  state  solution.  Water  is  often  incompressible  in 
applications  such  as  this.  With  assumption  1,  constant  average  speed  remains  valid  as  long 
as  the  pipe  diameter  is  constant.  Heat  transfer  being  the  only  form  of  loss  is  true  since  there 
are  no  work-extracting  elements.  Darcy  Friction  Factors  were  assumed  in  all  equation 
calculations  as  it  is  the  most  commonly  used  parameter  and  the  derivations  of  the 
equations2  imply  its  use.  Lastly,  assumption  6  is  valid  since  the  hope  is  to  have  a  highly 
conductive  pipe  material  for  this  application. 

If  this  last  assumption  did  not  hold,  the  left  hand  side  of  (9)  would  have  an  extra  thermal 
resistance  term.  This  term  would  be  a  thermal  resistor  added  in  series  to  the  heat  transfer 


2  Mills,  A.  (1999).  Basic  Heat  and  Mass  Transfer  (2nd  ed.).  Los  Angeles,  CA:  Pentice  Hall. 
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coefficient  term  (which  is  representing  a  conductance).  This  would  change  the  equation  in 
the  following  way: 


Nu  k 


D 


=  > 


ln((P+2t)/P) 


Nu k  2nkpl 

where  t  is  the  thickness  of  the  pipe  and  kp  is  the  thermal  conductivity  of  the  pipe.  This 
makes  the  relationship  between  l,  D,  and  the  heat  loss  more  complicated: 

,  _  FA+f (F A)2  +4FBCDETT 
_  2CDEn  '  ' 

A  =  2nDkv,  B  =  Nu  /dn((D  +  2 t)/D),  C  =  2nkvNu  k 


E  =  (T1  +  T2-  2Ts),  F  = 


pvnD2 


cM  -  t2) 


The  one  major  thing  that  is  uncertain  is  how  the  soil  temperature  in  contact  with  the  pipe 
is  to  be  modeled.  The  following  example  assumes  a  constant  soil  temperature.  With  all  the 
assumptions  present,  the  values  found  with  this  analysis  are  low  end  estimates. 


Example: 

For  this  project,  it  is  expected  that  there  will  be  a  large  flow  rate  of  water.  To  estimate 
this  number,  a  baseline  of  22  *  10 e gal/ day  was  suggested.  This  in  SI  units  is: 

22  *  106  — 

day 

Volume  flow  rate  is  related  to  speed  and  diameter  by: 


1  day 

1  hr 

1  min 

0.00378541178  m3 

24  hr 

60  min 

60s 

1  gal 

V  =  Av  = 


nDz 


V  =>  V  = 

tcD2 

13.64 

s 


=  0.9638778  — 

s 


4V 


7TL 

0.964  — 

s 


Choosing  a  pipe  diameter  of  D  =  0.3m  requires  v 

The  remaining  data  needed  for  this  problem  are  calculated  at  room  temperature  (i.e. 

300 K)  for  salt  water  (25%  by  mass  NaCl  solution3)  since  the  temperature  changes  are  not 
expected  to  be  drastic  enough  for  property  changes  and  salt  water  is  the  possible  working 
fluid: 

p  =  1178^,  p  =  1.65  *  10“3  ^f,  cv  =  3310—,  k  =  0.470—,  Tt  =  308 K, 

^  m3  ^  m2  p  kgK  mK ’  1 

T2  =  300 K 

With  this  information,  the  Reynold’s  number,  thermal  diffusivity,  kinematic  viscosity,  and 
Prandtl  number  are  calculated  is  calculated: 


Re  =  PJ£  =  =  2.92144  ,  X06,  „  =  JL  =  »  «»  ’"2 


1.65*10-3 


pCp 


1178*3310  s 


1.2053816  * 


10“7  — 


_  p  _  1.65*10~3  m2 
~  p~  1178  5 

use  (5)  in  (10): 

pD2vcp(r1-r2) 


.  . _  1.400679*10  6 

1.400679  *  10  6 — ,Pr  = 


l  = 


1.2053816*10- 


1178*0.32*  13.64*33 10(8) 


2(0.023  Re°-8Pr°-33fc)(r1+T2-27's)  2(0.023*(2.92144*106)°-8*(11.620212)°-33*0.470)(8) 


11.6202122215 


662  m 


3  Mills,  A.  (1999).  Basic  Heat  and  Mass  Transfer  (2nd  ed.).  Los  Angeles,  CA:  Pentice  Hall. 
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Pressure  Analysis  for  Pipe  System 

Description:  A  level  pipe  system  (ps)  in  a  radiator  like  layout  is  proposed  as  a  pipe  layout  to 
carry  water  in  a  radiant  cooling  project  design.  The  system  will  be  created  from  a  specified 
number  of  standard  straight  and  bent  sections,  /?  and  y  respectively 
Given: 


D  -  diameter  of  the  pipe  system 
v  -  average  speed  of  fluid  in  the  ps 
l  -  length  of  straight  part  of  system 

R  -  radius  of  curvature  of  the  ps  bend 

p  -  the  density  of  water 
Determine:  The  pressure  loss  throughout  the  system 
Model:  1)  Steady  Incompressible  Flow 

2)  No  working  components  in  the  system 

3)  Constant  average  speed 

4)  Constant  elevation 

5)  No  heat  or  internal  energy  converts  to  mechanical  energy 
Basic  Equations:  Control  Volume  Conservation  of  Energy 

Q-w  =  YtfepdV  +  fA{u  +  Y  +  gz  +  B3pv'dA  (1) 

Energy  per  unit  Mass 

y2 

e  =  u  +  —  +  gz  (2) 


Conservation  of  Mechanical  Energy 

y  +  («2  ^  ^)  +  A  (gz)  =  -h 


(3) 


Head  Loss  in  Straight  Pipes 
Head  Loss  in  Pipe  Bends 


i  _  r  l  V 

hiyt  —  f 
M  J  D  2 


hm  =  K"~ 


(4) 


(5) 


Analysis:  by  applying  the  assumptions  1  and  2  to  (1)  and  substituting  in  (2)  one  obtains  (3): 

v2  pV\ 

u  +  —  +  gz  H - jpv  •  d.A 

2  m  J 

1&2  ,  „  /  „2  _  (2)  .  .  /  „2 


IJAllg,  ClOOLAllipULWllO  X  0.1  IVa  ^ 

Q-W  =  ±jePdV  +  jA(i 


Q  =  fA  (u  +  y  +  gz  +  jjj)  pv  ■  dA  ->  Q  =  jA  (u  +  y  +  gz  +  jjj)  pv  ■  dA 
Q  =  m  (u2  -  ux  +  gz2  -  gz1  +  ^  +  fA  pv-dA-  JA  pv  ■  dA 


(  p2  Pi\  v2  vt 

=  m  ^u2  -  Ui  +  gz2  -  gzx  +  —  -  —)  +  m(a2  —  -  ax  — ) 

pAz) 


p  p  _  _ 

5  Q  r  N  7  /Ap  v22  v±2 

^--(U2-Ui)  =  -A=(7  +  a2  — -a1T 

The  derivation  was  shown  to  explain  the  a  parameters  as  constants  determined  by 
the  flow  type  that  allow  the  use  of  the  average  velocity  to  calculate  the  required  integrals 
of  the  fluid  flow.  This  expression  ultimately  says  that  the  changes  in  mechanical  energy 
of  a  fluid  flow  result  in  head  loss  when  there  is  heat  flow  from  the  system  or  energy 
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storage  in  the  fluid.  All  this  leads  to  using  (3)  as  the  base  equation.  Now  apply 
assumptions  3  and  4  to  simplify: 


—h  =  (^j  =>  A p  =  —hp  (6) 


Since  the  interest  is  just  in  the  magnitude  of  the  pressure  change  the  negative  will  be 
ignored  from  here  on  out. 

The  head  loss  term  ( h )  consists  of  two  terms:  one  caused  by  friction  in  the  straight 
pipe  components  and  one  caused  by  losses  in  the  bend  components.  Assuming  that  the  ps 
is  made  from  a  combination  of  /?  straight  components  and  y  bent  components,  the  total 
head  loss  will  just  be  a  sum  of  all  the  component  head  losses: 


»-(r*  T  +  ^Ft)  « 

Here  the  loss  coefficient  K  is  a  function  of  the  nondimensional  term  -  and  can  be  found 

D 

in  appropriate  tables  for  180°  bends.  The  friction  coefficient  is  a  nondimensional  term 
dependent  on  the  Reynold’s  number  (Re  =  and  Surface  Roughness  coefficient  (|) 
found  using  a  Moody  Plot  or  an  empirical  fitting  equation  like  the  Colebrook  equation 
(- jj  =  —2.0  log  For  laminar  flows,  there  is  an  analytic  expression  of 

64 

/  =  —  which  is  captured  with  the  Moody  plot.  Lastly,  the  speeds  to  use  in  the  equation 


are  the  speeds  at  the  respective  positions  in  the  pipe  where  the  loss  is  occurring.  With 
assumption  3,  the  speed  is  a  constant  value  shared  everywhere. 

From  these  conditions,  combining  (6)  and  (7)  yields 

hp  =  p(yk  + fifty  f  (8) 

Discussion:  This  analysis  provides  a  simple  relationship  between  the  pressure  loss,  length, 
diameter,  and  radius  of  bend  curvatures  for  a  desired  design  speed.  The  cost  of  computing 
this  simple  result  will  be  in  finding  K  and  /  for  the  different  geometries  and  materials  used. 
The  assumptions  used  in  this  analysis  are  justified.  The  steady  flow  aligns  with  the  desire 
for  a  steady  state  solution.  Water  is  often  incompressible  in  applications  such  as  this.  The 
no  working  components  assumption  was  made  so  that  this  analysis  can  help  size  pumps  to 
counteract  the  pressure  drop.  With  assumption  1,  constant  average  speed  remains  valid  as 
long  as  the  pipe  diameter  is  constant.  Constant  elevation  results  from  the  fact  that  the  ps  is 
supposed  to  be  level.  Lastly,  assumption  5  is  valid  since  there  is  a  natural  tendency  for 
mechanical  work  to  convert  to  heat  loss  or  an  increase  in  internal  energy.  When  the 
opposite  does  occur,  the  amount  converted  to  mechanical  energy  is  often  insignificant.  If 
the  exact  system  is  not  as  depicted  in  the  figure  above,  the  total  head  loss  term  will  equal  to 
the  sum  of  (4)  using  the  total  length  of  straight  pipe  and  terms  from  (5)  for  each  different 
type  of  non-straight  pipe  with  corresponding  K  values. 


Example: 

For  this  project,  it  is  expected  that  there  will  be  a  large  flow  rate  of  water.  To  estimate 
this  number,  a  baseline  of  1000 gal /min  was  suggested.  This  in  SI  units  is: 


1000  — 


1  min 

0.00378541178  m3 

60s 

1  gal 

=  0.06309196  — 

S 


0.06309  — 

S 


Volume  flow  rate  is  related  to  speed  and  diameter  by: 
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7ZDZ 


4V 


V  =  Av  =  — v  =>  r?  =  — - 

4  7TD2 


Choosing  a  pipe  diameter  of  D  =  0.3 m  requires  v  «  0.8925 


m 
s  * 


The  length  of  the  Thompson  Drydock  is  about  259  m4.  Assuming  that  the  project  would 
not  span  an  entire  drydock  length,  the  straight  length  will  be  set  to  /  =  100  m.  Then,  for 
convenience,  R  =  4.5 D  =  1.35m.  There  will  be  5  bends  and  6  straight  segments.  With  the 
five  bends  having  the  given  radius,  this  ps  still  fits  into  the  width  of  the  referenced 
Thompson  Drydock.  Last  of  the  given  data  would  be  the  fluid  properties  of  water:  p  = 


kg 


Ns 


1000  ^  &  p  =  8.93  *  10"4 


m 


With  this  information,  the  Reynold’s  number  is  calculated: 

pvD  1000  *  0.892  *  0.3 


Re  = 


1-4 


3.00  *  105 


p  8.93  *  10" 

Since  the  diameter  of  the  pipe  is  quite  large  compared  to  the  possible  roughness  lengths  of 
the  pipe,  the  pipe  can  be  considered  smooth  (i.e.  -  «  0).  With  this  knowledge,  the  Moody 
chart  gives: 


/  =  0.0145 

5  R 

Using  the  loss  coefficient  chart  ,  K  «  0.3  for  this  example  (thanks  to  the  choice  of  -  = 
4.5). 

The  rest  is  evaluation  using  (8): 

/  100\  0.89252 

A p  *  1000  ^5  *  0.3  +  6  *  0.0145  -— ) - - - Pa  *  12.15  kPa 


4  http://www.oracleireland.com/Ireland/Countvs/antrim/Shipvard/thompson-dock.htm 

5  http://www.thermopedia.com/content/577/?tid=104&sn=1422 


16 


Appendix  B  -  Python  Notebook  Documentation 

In  [  ] :  #  Import  needed  packages  as  abbreviations 
import  matplotlib.pyplot  as  pit 
import  numpy  as  np 


In  [  ]: 


# CHANGING  the  values  will  mess  up  some  of  the  labeling  formats 
# Define  important  values  of  the  fluid  bulk  (assumed  to  be  same  as  surface) 


rho  =  1178; 
mu  =  1 .65*10** (-3) ; 
cp  =  3310; 
k  =  0.470; 

T1  =  300+8*5. /9.; 

T2  =  300.0001; 

D_eng  =  np.arange(10,25,2) ; 
D  =  D_eng*0.0254; 


#. Density  of  fluid  (kg/m~3) 

# Dynamic  viscosity  of  fluid  (N*s/m~2) 

# Specific  heat  of  fluid  (J/kg/K) 

# Thermal  conductivity  of  fluid  (W/m/K) 
itlnlet  temperature  designed  to  be  8  deg 

#  F  above  ambient  (K) 
it Outlet  temperature  (K) 

# Diameters  in  inches  for  fomatting 
UDiameters  of  pipe  converted  from 

#  inches  to  meters 


UCalculate  the  length  of  D 
Dlength  =  len(D) ; 

e  =  0; 

Ts  =  300; 
kp  =  0.04; 


UChar act  eristic  Roughness  length 
UTemperature  of  surrounding  soil 
ttThermal  Conductivity  of  the  pipe  (W/m/K) 


HPipe  thickness  of  schedule  1^0  pipe  converted  from  inches  to  meters 
it  10  12  lj.  16  18  20  22  24  (in)  diameter  pipe 

Hcorresponding  to  the  thickness 

t  =  [0. 37, 0. 41, 0. 44, 0. 5, 0. 56, 0.59, 0.59,0. 69] *np. full (Dlength, 0.0254); 


ttSet  important  values  for  pressure  drop 

kCoef  =  0.3 

lPipe  =  239 

bends  =  5 

strts  =  6 

rBend  =  4.5*D 


ttHead  loss  coefficient 
ULength  of  parking  lot  (meters) 
UNumber  of  180  degree  bends  in  system 
UNumber  of  straight  parts 
HRadius  of  the  bends  (meters) 


itSet  information  about  pricing 
it  10  12  14  16  18  20 

insulPrice  =  [17.26,26.67,35.06,44.19,66.75,78.13, 
] /np . full (Dlength , 0 . 3048) ; 

condPrice  =  insulPrice* 103. 34/26. 67 


condPrice[0]  =  88.2/0.3048 

costPwr  =0.12 
pump_eff  =0.82 


22  24.  (in)  diameter  pipe 

U corresponding  to  the  thickness 
83.37,83.37 

#£  per  meter  converted  from 
it  per  foot  estimated 
#£  per  meter  based  on  how  one  price 
it  compared  to  plastic  of  same 

H  diameter 

ttAdjust  the  cost  of  the  first  item 
it  to  what  has  be  quoted 

#£  per  kWh 

it Efficiency  of  the  pump  purposely 
a  chosen  poor 


ttCalculate  the  thermal  diffusivity  and  kinematic  viscosity 
a  =  k/rho/cp  tt(m~2/s) 
nu  =  mu/rho  #(m^2/s) 


UDefine  vector  of  various  flow  rates  desired 

VdotG  =  np.arange(l,26.05,0. 1)  Million  Gallons  per  day 

ttConvert  the  flow  into  SI  units 

Vdot  =  VdotG*10**6*0. 00378541178/24/3600  UMeters  cubed  per  second 

HConvert  to  gpm 

VdotGPM  =  VdotG* 10**6/24/60 

HRecord  the  lengths  of  volume  flow  vectors  for  easy  use  later 
length  =  len(Vdot); 

In  [3] :  matplotlib  inline 

In  [4]  :  itPlot  with  respect  to  volume  flow  rates 
iiCreate  figures  with  different  handles 
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figl  =  pit .figure(l) ; 
fig2  =  pit . figure (2) ; 
fig3  =  pit .figure (3) ; 
fig4  =  pit .figure (4) ; 
fig5  =  pit . figure (5) ; 
fig6  =  pit .figure (6) ; 
fig7  =  pit . figure (7) ; 
fig8  =  pit .figure (8) ; 

ttCreate  axes  for  all  of  these 
axl  =  figl .add_subplot( 111) ; 
ax2  =  fig2.add_subplot(lll) ; 
ax3  =  fig3.add_subplot(lll) ; 
ax4  =  f ig4. add_subplot (111) ; 
ax5  =  fig5.add_subplot(lll) ; 
ax6  =  fig6.add_subplot(lll) ; 
ax7  =  fig7.add_subplot(lll) ; 
ax8  =  fig8.add_subplot(lll)  ; 

# Create  twin  axes  to  display  both  metric  and  english  units  side  by  side 

axil  =  axl.twinyO; 
ax22  =  ax2.twiny() ; 
ax33  =  ax3.twiny(); 
ax44  =  ax4.twiny(); 
cix55  =  ax5.twiny(); 
ax55  =  ax5.twiny(); 
ax66  =  ax6.twiny(); 
ax77  =  ax7.twiny(); 
ax88  =  ax8 . twiny () ; 

# Create  a  variable  for  consistent  formatting  later 

#The  when  variable  tells  what  power  cost  graphs  to  graph  by  setting  the  number  of 

#  initial  initial  vector  values  not  to  be  graphed.  Since  vectors  are  set  from 

#  small  to  large  diameter,  the  smaller  diameter  pipes  are  taken  out  first. 

when  =  0; 

#The  whenStop  variable  records  the  last  vector  position  that  the  pressure  graphs 

#  are  shown  to.  If  anything  other  than  -1  it  is  safest  to  take  of  the  formatting 

#  of  the  second  x  axis  designated  with  [*****] 

whenStop  =  -1; 

#. Looping  for  plotting 

for  num  in  xrange(0,Dlength) : 

# Create  the  velocities  for  the  graphs  (m/s)  from  flow  rate  and  diameter 
vels  =  4*Vdot/np.pi/D[num] **2; 

# Create  Reynolds  numbers  for  the  flows 
Re  =  rho*D[num] *vels/mu; 

# Create  the  Prandtl  Number  for  the  flows 
Pr  =  nu/a; 

#Check  if  Pr  is  above  0.5,  if  this  is  false  the  equations  used  are  invalid 
print  Pr  >  0.5 

#Preallocate  an  empty  vector  of  Nusselt  Numbers 
Nu  =  np. empty (length) ; 

^Calculate  the  nusselt  numbers 
for  dum  in  xrange (0, length) : 

# Determine  which  way  to  calculate  the  Nusselt  Number 
if  Re [dum]  <  3000: 

Nu [dum]  =  3.66; 

elif  Re [dum]  <  10**4:  #this  used  to  be  10**6  but  is  now  just  set  to  when 
#  the  other  equation  is  not  applicable 
# Approximate  Friction  factor  adds  2%  error 

fr  =  (-1.8*np.logl0((e/D[num]/3.7)**l.ll+6.9/Re[dum]))**(-2) ; 

Nu  [dum]  =  (  (Re  [dum]  -10**3)  *Pr*f  r/8)  /  (1+12 . 7*  (Pr**  (2/3)  -1)  *  (f  r/8)  **  (1/3)  )  ; 
else: 

Nu [dum]  =  0 . 023* (Re [dum] **0 . 8) * (Pr**0 . 33) 

# Calculate  the  length  1  required  to  get  the  temperature  drop 

11  =  np . divide (( (rho*D [num] **2*vels*cp)* (T1-T2) ) , ((2*Nu*k*(Tl+T2-2*Ts)))) ; 

# Calculate  length  2  using  chunks 


#For  conductive  length 
#For  insulative  len 
#For  Pressure  drop  in  conductive 
#For  Pressure  drop  in  insulative 
#For  Initial  Cost  conductive 
#For  Initial  Cost  insulative 
#For  Power  cost  conductive 
#For  Power  cost  insulative 
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A  =  2*np.pi*D[mim]  *kp; 

B  =  np.log((D[num]+2*t  [num]  )/D[num]  )*Nu*k; 

C  =  Nu*k*2*np.pi*kp; 

E  =  Tl+T2-2*Ts; 

F  =  (rho*D [num] **2*vels*cp) * (T1-T2) /2*np . pi ; 

12  =  (F*A+np .  sqrt  (  (F*A)  **2+4*E*C*np .  pi*D  [mm]  *F*B)  )  /  (2*E*C*np .  pi*D  [num]  )  ; 

#Preallocate  an  empty  vector  of  friction  factors  for  pressure  calculation 
f  =  np. empty ( length) ; 

# Calculate  the  friction  factors 
for  dum  in  xrange (0, length) : 

# Determine  which  way  to  calculate  the  friction  factor 
if  Re  [dim]  <  3000: 

f [dum]  =  64/Re [dum] ; 
else: 

# Approximate  Friction  factor  adds  2%  error 

f [dum]  =  (-1 .8*np . Iogl0( (e/D [num] /3.7)**1 . 11+6 . 9/Re [dum] ) )**(-2) ; 

# Calculate  the  pressure  drops  in  bar  //feet  if  the  comments  are  taken  avay 
pLossl  =  rho*(bends*kCoef+f*(ll-bends*np.pi*2*rBend[num] )/D[num] ) *vels**2/2/10**5\ 
#*10**5/rho/9 . 8/0. 30/8;  #What  to  uncomment  starts  if  gone 
pLoss2  =  rho*(bends*kCoef+f*(12-bends*np.pi*2*rBend[num] )/D[num] ) *vels**2/2/10**5\ 
#*10**5/rho/9. 8/0. 30/8;  #What  to  uncomment  starts  if  gone 

# Calculate  the  cost  for  the  different  materials  by  calculating  the  power  required 
pcostDayl  =  pLossl*10**2*Vdot/pump_eff ;  ttunits  kPa*m''3/s 
pcostDay2  =  pLoss2*10**2*Vdot/pump_eff ;  #units  kPa*m"3/s 

#Plot  the  11  data 
axl .plot (Vdot ,11) 

# Annotate  along  the  way 
pt  =  -1 

axl .  annotate (M,/,d  in"  °/,D_eng [num]  ,  xy=(Vdot  [pt]  ,  11  [pt]  +10)  .ha^  center  ’ ) 
pit .hold 

tiPlot  the  12  data 

ax2 . plot (Vdot , 12 , f igure=f ig2) 

#Label  the  graphs  with  certain  diameters 
if  num  ==  Dlength-1: 
pt  =  0 

ax2 . annotate ("%d"  %D_eng [num] ,  xy= (Vdot [pt] ,  12 [pt] +20) ,ha=  * right ’ ) 
pit .hold 

#Plot  the  pressure  data 
if  num  >=  when: 

ax3 .plot (Vdot [0 :whenStop] .pLossl [0 :whenSt op] ,f igure=f ig3) 
pit .hold 

ax4. plot (Vdot [0:whenStop] ,pLoss2 [0 :whenStop] ,f igure=f ig4) 
pit .hold 

# Annotate  the  pressure  data 
if  num  ==  0: 
pt  =  120 

ax3 .  annotate("%d  in"  °/,D_eng [num]  ,  xy=(Vdot  [pt]  ,  pLossl  [pt]  )  .ha^ right  ’ ) 
ax4. annotate ("®/.d  in"  °/,D_eng [num]  ,  xy=(Vdot  [pt]  ,  pLoss2  [pt]  +10)  ,ha=} right  ’ ) 
pit .hold 

#Plot  the  initial  cost  of  stuff  -  conductive 
intCostl  =  condPrice [num] *11 ; 
ax5. plot ( Vdot ,intCostl) ; 

#. Annotate  along  the  way 
pt  =  -1 

2lx5. anno t ate ("°/0d  in"  °/,D_eng [num]  ,  xy=(Vdot  [pt]  ,  intCostl  [pt] +6000)  ,ha=’ center  ’  ) 
pit .hold 

#Plot  the  initial  cost  of  stuff  -  insulative 
intCost2  =  insulPrice [mm] *12; 
ax6 . plot (Vdot , intCost2) ; 
it  Annotate  along  the  way 
if  num  ==  0: 
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pt  =  0 

ax6. anno tat e("°/.d"  °/.D_eng[num]  ,  xy=(Vdot  [pt]  ,  intCost2 [pt]  -3000)  ,ha=’ right  ’  ) 
pit .hold 

if  num  >=  when: 

#Plot  the  power  cost  of  conductive 
ax7 . plot (Vdot [0 : whenStop] , pcostDay 1 [0 : whenStop] ) ; 
pit .hold 
if  num  ==  0: 
pt  =  120 

ax7.annotate(",/,d  in"  °/«D_eng [num]  ,  xy=(Vdot  [pt]  ,  pcostDayl  [pt] +200)  , 
ha=’ right  * ) 

ax8. annotate ("°/,d  in"  °/0D_eng [num]  ,  xy=(Vdot  [pt]  ,  pcostDay2 [pt]  +200) , 
ha= ’ right ’ ) 

#Plot  the  power  cost  of  insulative 

ax8 . plot (Vdot [0 : whenStop] , pcostDay2 [0 : whenStop] ) ; 
pit .hold 

# Labels  for  the  graphs 

# Create  number  of  columns  to  make  lengend  not  over  lap  data 
ncol  =  2 

#[*****]  Label  the  always  the  same  x  axes 
for  num  in  xrange(0,8): 

funstr  =  "ax°/,d. set _xlabel(’ Volume  Flow  Rate  (nT3/s)’)"  °/,(num+l) 

funstr2  =  "ax°/,d°/0d.  set  _xlabel(’  Volume  Flow  Rate  (gpm)’)"  %(  (num+1)  ,  (num+1)) 

funstr3  =  "axyd'/.d.  set_xlim(VdotGPM[0]  ,VdotGPM  [-1]  )  "  °/,( (num+1)  ,  (num+1)) 

eval (funstr) 

eval(funstr2) 

eval(funstr3) 


# Conductive  Pipe  Required  Length  Figure 
axl .set_ylabel(’ Required  Pipe  Length  (m)’) 

axl .set_title(’ Length  vs.  Volume  Flow  Rate  in  Highly  Conductive  Pipe\n\n’) 

# Insulative  Pipe  Required  Length  Figure 
ax2.set_ylabel(’ Required  Pipe  Length  (m)’) 

ax2 . set .title ( ’Length  vs.  Volume  Flow  Rate  in  Insulating  Pipe\n\n’) 
ax2.1egend(D_eng[np.arange(when,Dlength,l)] ,  loc  =  0,  ncol=ncol, 
title= ’Diameter  (in)’) 

# Conductive  Pipe  Pressure  Loss  Figure 
ax3.set_ylabel(’ Pressure  Loss  (bar)’) 
ax3.set_title(’ Pressure  Loss  in  Conductive  Pipes\n\n’) 
ax3. legend (D_eng[np.arange(when,Dlength,l)] ,  loc  =  0,  ncol=ncol, 
title=,Diameter  (in)’) 

# Insulative  Pipe  Pressure  Loss  Figure 
ax4.set_ylabel( ’Pressure  Loss  (bar)’) 
ax4 . set_title( ’Pressure  Loss  in  Insulative  Pipes\n\n’) 
ax4. legend (D_eng [np . arange (when, Dlength, 1)] ,  loc  =  0,  ncol=ncol, 
title=’ Diameter  (in)’) 

# Conductive  Pipe  Initial  Costs  Figure 
ax5 . set_ylabel( ’Cost  ($)’) 
ax5.set_title(’Pipe  Cost  -  Conduct ive\n\n ’ ) 

# Insulative  Pipe  Initial  Costs  Figure 

ax6.set_ylabel(’Cost  ($)’) 

ax6 . set_title( ’Pipe  Cost  -  Insulative\n\n’ ) 

# Conductive  Pipe  Power  Cost  Figure 
ax7 . set_ylabel( ’Power  (kW)’) 
ax7.set_title(’ Power  Cost  -  Conduct ive\n\n’ ) 

clx7.  legend  (D_eng[np.  arange  (when,  Dlength,  1)  ]  ,  loc  =  0,title= ’Diameter  (in)’, 
ncol=ncol) 

# Insulative  Pipe  Power  Cost  Figure 
ax8.set_xlabel(’ Volume  Flow  Rate  (nT3/s)’) 
ax8 . set_ylabel( ’Power  (kW)’) 
ax8.set_title(’ Power  Cost  -  Insulative\n\n’ ) 

ax8. legend (D_eng [np . sircinge (when, Dlength,  1) ]  ,  loc  =  0,title= ’Diameter  (in)’, 
ncol=ncol) 
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#show  the  plots 
pit .  showO 

In  [14] :  #Plot  with  power  cost  and  initial  cost  with  respect  to  each  other 

# Create  the  design  variable  for  the  total  flow  rate  to  be  supplied 
totVdot  =  26*10**6*0.00378541178/24/3600  #Twenty-six  million  gpd 

#  converted  to  m^S/s 


# Create  figures  with  different  handles 

fig9  =  pit .figure (9) ;  #For  Power  cost  vs  initial 

#  cost  conductive 


#Create  axes  for  all  of  these 
ax9  =  fig9.add_subplot(lll) ; 

#Create  Color  Code 

color  =  [’b’ , 'g' , ’k’ , ’m* , ’r* , ’c’] 

# Establish  the  array  to  create  the  number  of  modules 
m  =  2**np.arange(0,5) 

# Initialization  step  for  the  vectors 

pcostl  =  np . empty (Dlength) ; 
pcost2  =  np. empty (Dlength) ; 
intCostl  =  np . empty (Dlength) ; 
intCost2  =  np . empty (Dlength) ; 

# Looping  for  plotting  conductive 
for  dum  in  xrange(0,len(m) ) : 

for  mim  in  xrange(0, Dlength) :  #Loop  through  the  different  diameters 

# Calculate  the  velocity  for  the  individual  modulesgraphs  (m/s)  from 

#flow  rate  and  diameter 

vels  =  4*totVdot/np.pi/D[num] **2/m[dum] ; 

ttCreate  Reynolds  numbers  for  the  flows 
Re  =  rho*D[num] *vels/mu; 

#Create  the  Prandtl  Number  for  the  flows 
Pr  =  nu/a; 

# Calculate  the  Nusselt  Numbers 
Nu  =  0 . 023* (Re**0 . 8) * (Pr**0 . 33) 

# Calculate  the  length  1  required  to  get  the  temperature  drop 

11  =  np . divide (( (rho*D [num] **2*vels*cp)*(Tl-T2)) , ((2*Nu*k*(Tl+T2-2*Ts)))) ; 

# Calculate  the  friction  factors 

f  =  (-1 .8*np.loglO((e/D [num] /3.7)**1 . 11+6 .9/Re) )** (-2) ; 

# Calculate  the  pressure  drops  in  bar  //feet  if  the  comments  are  taken  avay 
pLossl  =  rho*(bends*kCoef+f*(ll-bends*np.pi*2*rBend[num] )/D[num] )*vels**2/2/10**5\ 
#*10**5/rho/9 .8/0 .30/8;  #What  to  uncomment  starts  if  gone 

# Calculate  the  cost  for  the  different  materials  by  calculating  the  power  required 
pcostl  [num]  =  pLossl*  10**2*totVdot/pump_eff ;  ttunits  kPa*m'‘3/s 

^Calculate  the  cost  of  the  pipe  in  the  required  length 
# Conductive 

intCostl [num]  =  condPrice [num] *ll*m[dum] ; 

#Plot  the  points 

ax9.  loglog(intCostl , pcostl ,  ’‘/s-s’  ‘/color  [dum]  )  ; 

# Looping  for  plotting  insulative 
for  dum  in  xrange(0,len(m) ) : 

for  num  in  xrange(0, Dlength) :  #Loop  through  the  different  diameters 

# Calculate  the  velocity  for  the  individual  modulesgraphs  (m/s)  from  flow  rate  and 
#  diameter 

vels  =  4*totVdot/np.pi/D [num] **2/m[dum] ; 

# Create  Reynolds  numbers  for  the  flows 
Re  =  rho*D[num] *vels/mu; 
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# Create  the  Prandtl  Number  for  the  flows 
Pr  =  nu/a; 

# Calculate  the  Nusselt  Numbers 
Nu  =  0 . 023* (Re**0 . 8) * (Pr**0 . 33) 

# Calculate  length  2  (i.e.  for  insulative  pipes)  using  chunks 
A  =  2*np.pi*D[num] *kp; 

B  =  np . log( (D [num] +2*t [num] )/D[num] )*Nu*k; 

C  =  Nu*k*2*np.pi*kp; 

E  =  Tl+T2-2*Ts; 

F  =  (rho*D [num] **2*vels*cp) * (T1-T2) /2*np . pi ; 

12  =  (F*A+np.sqrt((F*A)**2+4*E*C*np.pi*D[num] *F*B))/(2*E*C*np.pi*D[num] ) ; 

# Calculate  the  friction  factors 

f  =  (-1 . 8*np . Iogl0( (e/D [num] /3 . 7) **1 . 11+6 . 9/Re) ) ** (-2) ; 

# Calculate  the  pressure  drops  in  bar  //feet  if  the  comments  are  taken  avay 
pLoss2  =  rho* (bends*kCoef +f * (12-bends*np . pi*2*rBend [num] ) /D [num] ) *vels**2/2/10**5 
#*10**5/rho/9. 8/0. 30/8;  #What  to  uncomment  starts  if  gone 
# Calculate  the  cost  for  the  different  materials  by  calculating  the  power  required 
pcost2[num]  =  pLoss2*10**2*totVdot/pump_eff ;  #units  kPa*m'‘3/s 

# Calculate  the  cost  of  the  pipe  in  the  required  length 
# Insulative 

intCost2 [num]  =  insulPrice [num] *12*m[dum] ; 

#Plot  the  points 

ax9.1oglog(intCost2,pcost2,  ’’/.s-'* ’  ’/.color  [dum]  ) ; 

# Labels  for  the  graphs 
#Label  the  always-the-same  axes 
ax9 . set_xlabel( ’ Initial  Cost  ($)’) 
ax9.set_ylabel(’ Power  Cost  (kW)  ’  ) 

#Show  the  legend 

ax9. legend (m,loc=0,title="  Number  of  Modules\n  Triangle  =  Insulative\ 

\n  Square  =  Conductive") #, ncol=ncol) 


#Add  title 

ax9 . set .title ( ’Pipe  Cost  Comparison’) 

#Label  a  arrow  to  show  direction  of  increasing  D 

ax9 . annotate ("Increasing  D",  xy=(10**6.5,10**2) ,  xytext=(10**5 .7, 10**3. 7) , 
arrowprops=dict (arrowstyle= ’->’), ha= ’ center ’ ) 


pit . show() 


In  [6]  :  # Create  A  plot  without  purely  bad  options 

#Plot  with  respect  to  power  cost  and  initial  cost 
# Create  the  design  variable  for  the  total  flow  rate  to  be  supplied 

totVdot  =  20*10**6*0.00378541178/24/3600  # Twenty  million  gpd  converted  to  m~3/s 

# Create  figures  with  different  handles 

figlO  =  pit .figure (10) ;  #For  Power  cost  vs  initial  cost 

#  conductive 


# Create  axes  for  all  of  these 
axlO  =  figlO. add_ subplot (111) ; 

# Create  Color  Code 

color  =  [’b’ , ’g’ , ’k’ , ’m’ , ’r’ , ’c’] 

# Establish  the  array  to  create  the  number  of  modules 
m  =  2**np.arange(0,5) 
base  =  np. arange(0,5) 

ttCreate  a  5  by  len(m)*Dlength  vector  where  first  column  is  initial  pipe  cost,  second 

#  column  is  power  cost,  third  is  a  material  identifier  (0=conductive, 

#  l=insulative) ,  fourth  is  diameter  of  the  pipe  used,  and  fifth  is  the  number  of 

#  modules  used  (m) 

# Initialization  step  for  the  vector 
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ptlist  =  np. empty ( [5,len(m)*Dlength*2] ) ; 

# Create  a  counter  that  starts  at  0 
counter  =  -1; 

# Looping  to  fill  the  vector  with  points 
for  dum  in  xrange(0,len(m)) : 

for  num  in  xrange(0,Dlength) :  #Loop  through  the  different  diameters 

# Calculate  the  velocity  for  the  individual  modulesgraphs  (m/s)  from  flow  rate 
#  and  diameter 

vels  =  4*totVdot/np.pi/D[num] **2/m[dum] ; 

#Create  Reynolds  numbers  for  the  flows 
Re  =  rho*D[num] *vels/mu; 

# Create  the  Prandtl  Number  for  the  flows 
Pr  =  nu/a; 

# Calculate  the  Nusselt  Numbers 
Nu  =  0 . 023* (Re**0 . 8) * (Pr**0 . 33) 


# Calculate  the  length  1  required  to  get  the  temperature  drop 

11  =  np . divide ( ( (rho*D [num] **2*vels*cp) * (T1-T2) ) , ( (2*Nu*k* (Tl+T2-2*Ts) ) ) ) ; 

# Calculate  length  2  (i.e.  for  insulative  pipes)  using  chunks 
A  =  2*np.pi*D[num] *kp; 

B  =  np.log((D [num]+2*t [num] )/D[num] )*Nu*k; 

C  =  Nu*k*2*np . pi*kp ; 

E  =  Tl+T2-2*Ts; 

F  =  (rho*D [num] **2*vels*cp) * (T1-T2) /2*np . pi ; 

12  =  (F*A+np . sqrt ( (F*A) **2+4*E*C*np . pi*D [num] *F*B) ) / (2*E*C*np . pi*D [num] ) ; 

# Calculate  the  friction  factors 

f  =  (-1 . 8*np . Iogl0( (e/D [num] /3.7)**1 . 11+6 . 9/Re)) **(-2) ; 

# Calculate  the  pressure  drops  in  bar  //feet  if  the  comments  are  taken  avay 
pLossl  =  rho*(bends*kCoef +f*(ll-bends*np.pi*2*rBend[num] )/D[num] ) *vels**2/2/10**5\ 
#*10**5/rho/9 . 8/0 . 30/8;  #What  to  uncomment  starts  if  gone 
pLoss2  =  rho*(bends*kCoef+f*(12-bends*np.pi*2*rBend[num] )/D[num] ) *vels**2/2/10**5\ 
#*10**5/rho/9 .8/0 .30/8;  #What  to  uncomment  starts  if  gone 

# Increment  the  counter 
counter  =  count er+1; 

#print  counter 


# Calculate  the  costs  for  conductive 

ptlist [0 , counter]  =  condPrice [num] *11 *m [dum] ; 

ptlist  [1,  counter]  =  pLossl*10**2*totVdot/pump_ef  f ;  #units  kPa*m''3/s 

ptlist [2, counter]  =  0; 

ptlist [3 , counter]  =  D_eng [num] ; 

ptlist [4, counter]  =  m[dum] ; 

# Increment  the  counter 
counter  =  count er+1; 


# Calculate  the  costs  for  insulative 

ptlist [0, counter]  =  insulPrice [num] *12*m[dum] ; 

ptlist [1 , counter]  =  pLoss2*10**2*totVdot/pump_ef f ; 

ptlist [2, counter]  =  1; 

ptlist  [3,  counter]  =  D_eng[num]  ; 

ptlist  [4, counter]  =  m[dum]; 


#units  kPa*m''3/s 


#Get  rid  of  for  sure  bad  options,  those  that  are  to  the  top  right  of  other  points 

# Create  a  temp  variable 

temp  =  ptlist; 

print  np. shape (ptlist) 

for  num  in  xrange(0,len(m)*Dlength*2) : 

for  dum  in  xrange(0,len(m)*Dlength*2) : 

if  temp[0,num]>=ptlist[0,dum]  and  temp  [1 , num]  >=ptlist  [1  ,dum]  : 

#Make  the  (joint  equal  to  the  smaller  point  if  it  is  in  the  top  right 
#  relative  to  that  point 
temp[: , num]  =pt list  [:  ,dum]  ; 


23 


#Find  the  unique  points  in  temp  and  their  indices  to  label  the  graph  appropriately 
nothing,  indeces  =  np. unique ( temp [0, :] , return. index=True) 

#Plot  the  remaining  points  with  annotations 
# Create  an  evenness  counter 
even  =  0; 

for  num  in  indeces: 
even  =  even  +  1 ; 

#Choose  the  symbol  to  be  used 
if  temp[2,num]==0: 

symbolUsed  =  ’'/.ss’  %  (color  [np.  select  (temp  [4,  num]  ==m,  base)]); 
else : 

symbolUsed  =  ’‘/.s'*’  '/.(color  [np.  select  (temp  [4,  num]  ==m,  base)]); 

#Plot  the  point 

axlO . loglog (temp  [0, num] , temp [1, num] .symbolUsed); 

#. Annotate  the  point  with  diameter 
if  even'/,2  ==  1 : 

axlO.  annotate ("°/.d"  '/.temp [3, num]  ,  xy=  (t emp  [0, num]  , temp [1, num]  )  , 
xytext=(temp  [0,num] *2, temp [1 ,num] *1 . 75) , 
arrowprops=dict (arrowstyle= ’ - ’ ) , ha= ’ center ’ ) 

else : 

axlO .  annotate ("'/.d"  '/.temp [3, num]  ,  xy=  (temp  [0 , num]  , temp [1 , num]  )  , 
xytext= (temp [0 , num] / 1 . 25 , temp [1 , num] /3) , 
arrowprops=dict (arrowstyle=  ’  - ’ ) ,ha= ’ center 1 ) 

#The  Legend  is  the  same  as  before  just  I  do  not  know  how  to  get  it  to  show 
#Label  the  graph 

axlO . set_title( "Best  Options  Available") 
axlO. set_xlabel(’ Initial  Pipe  Cost  ($)’) 
axlO . set_ylabel( ’Power  Cost  (kW)’) 


24 


Appendix  C  -  Relevant  but  Less  Important  Curves  Generated 

Pressure  Loss  in  Conductive  Pipes 


Volume  Flow  Rate  (gpm) 

5000  10000  15000 


Fig.  Cl  -  Pressure  Loss  in  Conductive  Pipes 
Pressure  Loss  in  Insulative  Pipes 


Volume  Flow  Rate  (gpm) 

5000  10000  15000 


25 


Cost  ($)  Cost  ($) 


Pipe  Cost  -  Conductive 
Volume  Flow  Rate  (gpm) 


5000  10000  15000 


Volume  Flow  Rate  (m~3/s) 


Fig.  C3  -  Cost  of  Conductive  Pipe 
Pipe  Cost  -  Insulative 
Volume  Flow  Rate  (gpm) 


5000  10000  15000 


Volume  Flow  Rate  (m~3/s) 


Fig.  C4  -  Cost  of  Insulative  Pipe 
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Power  Cost  -  Conductive 
Volume  Flow  Rate  (gpm) 

5000  10000  15000 


Fig.  C5  -  Power  Cost  for  Conductive  Pipes 
Power  Cost  -  Insulative 
Volume  Flow  Rate  (gpm) 


5000  10000  15000 
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Power  Cost  (kW) 


Pipe  Cost  Comparison 
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Appendix  D  -  Price  Summary 


Diameters 

Insulative  Pipe  Cost* 

Conductive  Pipe  Cost** 

(in) 

($/ft) 

($/ft) 

10 

17.26 

88.20 

12 

26.67 

103.34 

14 

35.06 

135.85 

16 

44.19 

171.23 

18 

66.75 

258.64 

20 

78.13 

302.74 

24 

112.31 

435.17 

Notes  *  -  Prices  based  on  PVC  Piping  as  of  July  2015 

**  -  Non-hyperlinked  prices  based  on  ratio  b/w  the 
12in  diameter  prices 


Appendix  E  -  Estimated  Install  Costs  for  One  System  Size 


Pipe  24"  diameter  x  2276  Linear  ft 

$15.85 

Million 

Pumps  1.3  HP  1200  GPM  (16  at  $3,300  each) 

$0.05 

Million 

Excavation  of  60000  cubic  yrds 

$8.68 

Million 

10%  Redundancy  Costs 

$2.46 

Million 

$27.04  Million 


1  -  Price  Quoted  from  eBay  :  Molds  &  Plastic  Machinery,  Inc.  (2015). 

Aurora  Pump  -  Centrifugal  Pump,  Type  411  BF,  1200  GPM,  185 
Head  Feet  1750  RPM.  Retrieved  August  6,  2015. 

2  -  Price  Quoted  from  [5] 
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Symbol  List 

A  -  area  (for  system  and  pipe  system  cross  section) 

D  -  diameter  of  the  pipe  system  (ps) 

K  -  pressure  loss  coefficient  in  180°  bends 

P  -  power 

Q  -  heat  transfer  rate 

R  -  radius  of  curvature  of  the  ps  bend 

Re  -  Reynold’s  number 

Tave  -  average  fluid  temperature 

T±  -  inlet  design  temperature 

T2  -  outlet  design  temperature 

Ts  -  temperature  of  the  soil  surrounding  the  pipe 

V  -  volume  flow  rate 

W  -  work  rate 

e  -  system  energy  per  unit  mass 

/  -  Darcy  friction  factor 

g  -  acceleration  of  gravity 

h  -  total  head  loss  (energy  per  unit  mass) 

hM  -  major  head  losses 

hm  -  minor  head  losses 

l  -  length  of  straight  part  of  system 
k  -  thermal  conductivity  of  working  fluid 
kp  -  thermal  conductivity  of  an  insulative  pipe 
m  -  mass  flow  rate 

v  -  fluid  speed  in  the  ps 

v  -  average  speed  of  fluid  in  the  ps 

t  -  insulative  pipe  thickness 

z  -  elevation  relative  to  a  referance 
a  -  thermal  diffusivity  of  working  fluid 
/?  -  number  of  bends  in  the  ps 

A p  -  pressure  drop 

i]  -  pump  efficiency 

e  -  surface  roughness 

Y  -  number  of  straight  sections  in  the  ps 

g  -  dynamic  viscosity  of  working  fluid 

v  -  kinematic  viscosity 

p  -  the  density  of  working  fluid 


